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^ . We present a technique for derivation of a priori bounds for Gevrey-Sobolev norms of 

Ps| ! space-periodic three-dimensional solutions to evolutionary partial differential equations 
^ \ of hydrodynamic type. It involves a transformation of the flow velocity in the Fourier 
space, which introduces a feedback between the index of the norm and the norm of 
the transformed solution, and results in emergence of a mildly dissipative term. We 
illustrate the technique using it to derive finite-time bounds for Gevrey-Sobolev norms 
of solutions to the Euler and inviscid Burgers equations, and global in time bounds for 
the Voigt-type regularisations of the Euler and Navier-Stokes equation (assuming that 
^ ' the respective norm of the initial condition is bounded). The boundedness of the norms 
• ' implies analyticity of the solutions in space. 

1 Introduction 

^ I We suggest a new approach to derivation of bounds for analytic Gevrey-Sobolev norms 
■ (see their definition (6) below) of solutions to evolutionary partial differential equations 
. of hydrodynamic type, which exploits a feedback between the norm of the suitably trans- 
^ ! formed solution and the index of the norm. 

Whether three-dimensional solutions to the Euler equation remain smooth forever is 
O I an open question. If the initial flow velocity belongs to the Holder class C^~^°' (a > 
is arbitrary), the solution remains in it at least for a finite time [4]. Moreover, if ini- 
tially the flow is analytic in space, the solution is guaranteed to be analytic in a finite 
spatio-temporal region [1, 2, 3, 5, 6]; the loss of analyticity does not occur before the 
^ i solution ceases to be continuously differentiable [2]. These demonstrations involve con- 
I struction of bounds for iterated approximations of vorticity followed along characteristics, 
the knowledge of the Green function in the integral representation of the flow in the terms 
of vorticity, and application of methods of the theory of functions of complex variables. 
It is yet unknown, whether a singularity can develop at finite times. 

Foias and Temam [11] proved that viscosity can enhance the smoothness of space- 
periodic solutions to the Navier-Stokes equation (see also [10]): If initially the flow velocity 
belongs to the Sobolev space Hi(T^) (i.e., the space of space-periodic functions, whose 
derivatives of order 1 are square-integrable), then at any small time t > the flow has 
a bounded Gevrey-Sobolev norm. Their method was used subsequently to demonstrate 
that a similar result holds for space-periodic solutions to various modifications of the 
Navier-Stokes equation [16, 31]. The method relies on the presence of the viscous term 
in the equation, and hence is not directly applicable to the Euler equation. 
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We overcome the difficulty by performing a suitable transformation of the Fourier 
coefficients of the flow, which introduces a dependence of the index of a Gevrey-Sobolev 
norm on the norm of the transformed solution. The transformation results in emergence 
of a weakly dissipative operator in the modified equation and it becomes possible to 
proceed following the method of [11] (time dependence of the index giving rise to a 
weakly dissipative operator was employed in [22]). In Section 2 we use this technique 
to give a simple proof of the boundedness of a Gevrey-Sobolcv norm of a space-periodic 
three-dimensional solution to the Euler equation on a finite time interval. In section 3 
the same technique is applied to the Burgers equation. 

The three-dimensional Navier-Stokes equation can be regularised by introduction of 
suitable nonhnear [18, 19, 20] or linear terms, such as hyperviscosity [23, 24], or the time 
derivative of the Laplacian of the flow (proposed by O.A. Ladyzhenskaya at the Interna- 
tional mathematical congress in 1966) which can be regarded as a non-standard "dynamic" 
viscosity. Boundary value problems for the regularised equation of this type, known as 
the Navier-Stokes- Voigt equation and the Oskolkov equation in the Russian literature, 
were studied in [25, 26, 27, 8]. They describe the motion of non-Newtonian visco-elastic 
fluids with memory (characterised by an exponential decay of the rate of deformation 
under a constant stress) arising for a class of linear integro-diffcrcntial rheology equations 
[28, 29, 30]. While in the Voigt regularisation of the Navier-Stokes equation dynamic vis- 
cosity is positive, its negative values may be also physically sound [27]. Regularisation by 
progressive damping of high-wavenumber harmonics in some instances of occurrence of the 
flow velocity in the equation gives rise to the so-called Camassa-Holm or Navier-Stokes-a 
equation [12, 13, 16] supposedly describing Lagrangian-averagcd flows [15]; however, an 
additional term must be introduced for this interpretation to be accurate [32, 33]. 

Bounds for Gevrey-Sobolev norms were derived in [21] for solutions to the so-called 
Voigt regularisation of the Euler equation (this type of regularisation was suggested for 
the Euler equation in [8]) by the method, initially applied in [22] to the "lake equation" 
(i.e., a modifled two-dimensional Euler equation). We consider a Voigt-type regularisation 
of the three-dimensional Euler equation in Section 4 and show that a milder regularisa- 
tion suffices to estabhsh the boundedness of Gevrey-Sobolev norms of solutions. The 
boundedness of the norms is demonstrated for solutions to the Voigt-type regularisation 
of the Navicr Stokes equation in Section 5; the presence of diffusion enables us to further 
weaken the regularising term. Constructions of Sections 4 and 5 illustrate another aspect 
of our approach: dependence of the index of a Gevrey-Sobolev norm on the norm of 
the transformed solution allows to decrease the order of nonlinearity in the energy-type 
inequality and thus to obtain global (in time) bounds for the norm. 

In this paper, we consider zero-mean (note that the means of solutions are conserved 
in time) space-periodic three-dimensional solutions, assuming for simplicity that the el- 
ementary periodicity cell is the cube = [0, 27r]^. It is not difficult to generalise our 
analysis to encompass the case of arbitrary periods along the Cartesian axes. All the 
equations mentioned above normally involve a force f . Following a long-established tra- 
dition, for the sake of simplicity we consider only the illustrative equations for f = 0. 
Reinstating the forcing in the analysis is also straightforward (provided the force has the 
right analyticity properties). 

Let us also remark on the accurate procedure for derivation of a priori bounds. The 
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bounds are derived below for Fourier-Galerkin truncations of the solutions. We do this 
by tacitly using the ODE's governing the evolution of the Fourier coefficients of the trun- 
cated solutions. The bounds involve norms of the initial conditions for the truncations. 
Assuming in their place (in general, larger) norms of non-truncated initial conditions, we 
obtain bounds for the truncated solutions, that are uniform in the number of harmonics 
retained in the Fourier-Galerkin truncation. By the standard arguments, these bounds 
remain valid for the norms of the non-truncated solutions. 



2 The three-dimensional Euler equation 

In this section we establish bounds for Gevrey-Sobolev norms of solutions to the Euler 
equation on a finite time interval [0,t^), assuming that the initial condition has a finite 
Gevrey-Sobolev norm. 

The motion of a perfect (inviscid) fiuid with velocity v under the action of a body 
force f is governed by the Euler equation 

dv 

_ + (v V)v = f + Vp. (1) 

Incompressibility of fluid implies the solenoidality of the flow: 

V • V = 0. (2) 

The initial (at t — 0) distribution of the flow velocity v^'"^) is prescribed. 
We expand a zero-mean space-periodic flow in into the Fourier series: 

V = E vne^" " (3) 

117^0 



(4) 



The flow is real, as long as 
and solenoidality (2) implies 

v„ • n = 0. (5) 

The norm || ■ ||q of a zero- mean space-periodic vector or scalar field v in the Sobolev 
space Hq{T^) is defined in the terms of its Fourier coefficients by the relation 

n^O 



where | ■ | denotes the Euclidean norm of a three-dimensional vector. For any cr > we 
define the Gevrey-Sobolev norm^ 

|w|^,,^E|wn|V-l"l|nr^. (6) 

117^0 



-"^ Functions, whose Gevrey-Sobolev norms are finite, are analytic (sec [22]), and hence the norm (6) 
might be called the analytic Gevrey-Sobolev norm. We could extend our analysis to non-analytic Gevrey 
classes of index a, < a < 1, whose norms involve the exponents exp(2cr|n|") instead of exp(2c7|n|) in 
(6). Under the transformation of equations that we employ a weakly dissipative operator emerges, whose 
symbol would then grow as |n|". Since the operator is the strongest for a = 1, we consider only this case. 
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Here the first index cr is a lower estimate for the radius of the region of analyticity of w 
around the real axis in the complex space. For a = 0, (6) defines the norm in H^{T^). 
Let I • \q denote the norm in the Lebesgue space Lq{T^). By the Sobolev embedding 
theorem [7, 35], for any positive q < 3/2 there exists a constant Cq such that for any 
function / e Hq{T^) 

l/l6/(3-2,) < CMU- (7) 

Let Vn denote the linear operator of projection of a three-dimensional vector on the 
plane normal to n 7^ 0: 

|np 

The evolution of Fourier coefficients of the flow is governed by equations 

^ + i X:(vk ■ (n - k)) P„v„_k = (8) 

obtained by substitution of the series (3) (and f = 0) into the Euler equation (1). 

Assuming | v^'"^) |o.^5_|_3/2 < 00 for some positive a and s < 1/2, we consider a transfor- 
mation 

Mt) = w„(i) exp (-/3|n|||w(x,i) 117^3/2) , (9) 

w(x,^) = $:wn(^)e■"•^ 

n 

where /3 and e < 2 are positive constants. The transformation involves solving the system 
of nonlinear equations (9) in Wn(t); the solution takes the form 

w„(i) = v„(t) exp(V'(t)|n|), 

where the quantity ■ip{t) > satisfies the equation 

It has a unique solution for any t > 0, because the l.h.s. is a continuous monotonically 
increasing function of ip (as discussed in the introduction, the a priori estimates are derived 

for Fourier-Galerkin truncations of the flow, whereby the sum (3) is assumed to involve 
a finite number of terms). If /3 < o"|v*^^'^^|^ s+3/2' then ||w(x, t)||s+3/2 is bounded at t = 
(uniformly over the number of terms in the truncated sum (3)). Our goal is to show that 
w is bounded in Hs_^.3/2{T^) on some finite-length time interval [0,t*). 
Substitution (9) transforms (8) into 

^ + ^£|n|||w(x,t)|| 7^3/2 w„ ^ 1 1 w (x, i ) 1 1 ,+3/2 
= -iE(wk • (n - k)) Kw„_kexp (/3||w(x,i) 11-3/2 (|n| - |k| - |n - k|)) . (10) 

k 

Scalar multiplying this equation by Inp"*"^* Wn and considering the real part of the sum 
over n, we find 
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|(^l|w(x,t)||?+, + ^||w(x,t)||^-,j 
= Im^(wk • (n - k))(w„_k -^11112+2^ exp {(3\\w{^,t)\\;+3/2 M - |k| - |n - k|)) . 

k,n 

(11) 

We call (11) the energy balance equation. In construction of estimates for the sum aris- 
ing from the nonlinear term we follow, with some variations, the approach [9]. By the 
inequality for sides of a triangle, the exponential function in the r.h.s. of (11) docs not 
exceed 1. We will show now that the absohitc value of the r.h.s. of (11) is bounded by 
(D5/2)||w||^_,_3^2) where the constant Dg is independent of w {Dg ^ oo for s ^ 0). 
By virtue of the inequality 

|n|^+V2< |n-k|^+V2 + |k|^+V2, 
valid for < s < 1/2, the r.h.s. of (11) does not exceed 

Y: (|wk||n - k|^+=^/2|w„_k| + |k|^+V2|^^||n _ k||w„_k|) |n|^+^/2|w„| 

n,k 

= {2nr' 1^^^ (/o(x)/,+3/2(x) + /,+i/2(x)/i(x)) /.+3/2(-x) dx, (12) 
where scalar functions fq are defined as the Fourier series 

/,(x,i) = ^|w„(i)||n|V"-. (13) 

n 

By the Cauchy-Buniakowski-Schwarz inequality, for any / = J2n fn^^'^'^: 

\ 1/2 



I/I 



<Y.\^r"^"\fM~^''-'^"^ < 11/11^+3/2 (ei^I"'"'^) 



The second factor in the l.h.s. of this inequality, which we denote by c^, is finite for any 
s > (it tends to infinity, when s — > 0). Consequently, using Holder's inequality, the 
Sobolev embedding theorem (see (7)) and Parseval's identity, we find a bound for the 
r.h.s. of (12): 

(27r)~^ (|/5+3/2|2max|/o| + |/,+i/2y/i|3) I/.+3/2I2 

< (27r)~^ (c5||/o||^+3/2 + C'lCl/2||/o||5+3/2||/o||3/2) II/0II5+3/2 < ^||w||^+3/2 

for Dg = {cg + CiCi/2)/(47r3). Thus, we obtain from (11) 

|(l|w(x,t)||?+, + A||w(x,t)||^-2) <^.l|w(x,t)||^+3/2, 

where 

A = 2l3e/{2-e). (14) 



5 



Hence, for ^ = ||w(x, + A||w(x, i)!!^-/^, 

^<D, ^-3/(2-^) £3/(2-e) - -r ^ < L'.M-^/^^-.) 

where 9 = {1 + e)/{2 — e). Consequently, the bound 

^ < {{^\t=or - DM-'^^'-'hy'^' (15) 

holds for 

t<U = {D,e)-'A^/(-^-'\^\t=o)-'. (16) 

For such t, 

||w(x,i)||,+3/2 < V{t) ^ {{A/^\t=of - DM/aY"^"'''^ . (17) 

Assuming in 

e|t=o= ||w(x,0)||?+, + A||w(x,0)||^-/2 (18) 
the norms of the non-truncated vector field w(x, 0), we obtain bounds for ||w(x, t)||s+3/2 
and that are uniform in the number of harmonics in the Fourier-Galerkin truncation 
of a solution. Inequality (17) and transformation (9) imply 

|v(x,t)|^^-.,3+3/2 < \^{^.t)\f)\\^{^,t)\\:U/2,s+i/2 = l|w(x,t)||,+3/2 < if{t). 

We have therefore proved 

Theorem 1. Let initial condition v*^™) of a solution to the force-free Euler equation 
have a finite Gevrey-Sobolev norm |v("^) s+3/2, where < s < 1/2. For < t < the 
solution satisfies the bound 

|v(x,t)|^^(t)-.,,+3/2 < v?(t), (19) 

where < /3 < o'|v*^"^''|^ s+3/2' < e < 2, and Lp{t) are defined by relations (16)-(18), 
and w(x, 0) is the result of application of transformation (9) to the initial condition v*^™). 

3 The inviscid Burgers equation 

In this section we show that our technique gives the same bound for the Gevrey-Sobolev 
norms of solutions to the inviscid Burgers equation 

^ + (v.V)v = f. 

The solenoidality of solutions is not required any more. For the sake of simplicity, we 
again consider only zero-mean space-periodic solutions for the force-free case f = 0, with 
the elementary periodicity cell being the cube T^. 

A solution is expanded into the Fourier series (3); the evolution of Fourier coefficients 
is now governed by equations 

^ + iX:(vk-(n-k))v„_k = 0. (20) 
After transformation (9) is applied, (20) becomes 
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-^+(5e\u\ ||w(x,t)||^^3/2W„— ||w(x,t)||,+3/2 
= -iE(wk • (n - k)) w„_kexp (^||w(x,i) 11^3/2 (|n| - |k| - |n - k|)) . (21) 

k 

Scalar multiplying this equation by Inp"*"^* Wn and considering the real part of the sum 
over n, we again obtain the energy balance equation (11). Since the orthogonality (5) was 
not used to derive from (11) the bound (17) yielding (19), the same derivation applies 
for the force-free inviscid Burgers equation, implying that the bounds (17) and (19) hold 
true for its solutions as well. In other words. Theorem 1 apphes hterally to solutions to 
the force-free inviscid Burgers equation. 



4 Voigt-type regular isat ion of the Euler equation 

In this section we present another illustration of our technique, demonstrating global (in 
time) boundedness of Gevrey-Sobolcv norms of solutions to the Voigt-type regularisation 
of the Euler equation. The name "Voigt regularisation" was proposed in [21] for s = 1 
(see (22) below). The boundedness was proved ibid, for s = 1 by the method [22]. 
Our technique seems simpler in application. We show that a milder regularisation, than 
considered in [21], suffices to guarantee global regularity of solutions. 

We consider three-dimensional solenoidal space-periodic zero-mean solutions (3) to the 
equation 

«'^(-V')'v + ^ + (v ■ V)v = f + Vp (22) 

for s > 0. For s = 1 this is the Voigt-regularised Euler equation, studied in [21]. For 
f = 0, in the Fourier space it takes the form 

(1 + «>r^)^ + i E(vk ■ (n - k)) P„Vn_k = 0. (23) 

We assume |v^™^ |o-_ <j+i/2 < oo for a u > and make the transformation 

v„(i) = w„(i)exp(-/3]n] ]||w(x,i)|||-^) , (24) 
w(x,^)=Ew„(^)e-^ 

n 

where ^ and e <2 are positive constants and 

ll|w|ir = E(l + «>r')|n||w„|2 (25) 

n 

defines a norm equivalent to || ■ ||s+i/2. Choosing 

< ^ < a(|vM|^,,/, + a^|v(^^)|^,,+,/,)^/\ (26) 
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we ensure that |||w(x, 0)||| is bounded uniformly over the number of Fourier harmonics 
preserved in Fouricr-Galcrkin truncations of the initial condition v^™\ Our goal is to 
derive bounds for |||w(x, 

Substitution (24) transforms (23) into 

(: + „>l^.)(^|.,,.|„lll|w||r-*4|||w|||) 

= -i E(^k • (n - k)) P„w„_kexp {/3\\\w\\\-' (|n| - |k| - |n - k|)) . 

k 

The real part of the sum over n of these equations scalar multiplied by Wn reduces to 

^(« + «^l|w||^ + A|||w|f-^) 
= 2Im^(wk • (n - k))(w„_k • ^) exp (/3|||w|||-^ (|n| - |k| - |n - k|)) 

k,n 

= Im53(wk • (n - k))(w„_k ■ ^) (exp (/3|||wr^ (|n| - |k| - |n - k|)) 

k,n 

-exp(/3|||w|||-^(|n-k|-|k|-|n|))) (27) 

(above, in one of the sums in k and n we have changed the index of summation n ^ k — n 
and rehed on the fact that w is real- valued and solenoidal, which stems from (5) and (4)). 

Note that |e" -e" | < | a' — a"\ for any negative a' and a". Consequently, the absolute 
value of the r.h.s. of (27) does not exceed 

^|wk||n-k||w„_k||wn|2/3|||w|||-^||n|-|n-k|| < ^'"^f ' / /o(-x)/i'(x) dx 

n,k ^'^ 

(here and in (27) we use the notation A and fg (13) introduced in Section 2) 

~^Tr^ ' ' 9/(2-30 I /l 1 18/(7+30 ^"^^f 

(by Holder's inequahty for a sufficiently small C > 0). By the Sobolev embedding theorem 

I /o 1 9/(2-30 < C'5/6+cll/o||5/6+0 |/l| 18/(7+30 - C'l/3-C/2 ||/l || 1/3-C/2 

(see (7)). By virtue of (13) and by Holder's inequahty, (28) is bounded by 

/3C5/6+cCl/3-C/2ll|w||r^ 2 ^ /5C'5/6+cCl/3-C/2lll'^lir^l ||1+6C || ||2-6C 

l|w||5/6+cl|w||4/3-C/2 < l|w||5/6+cl|w||4/3+^. 

Hence, for s > 5/6 + C and s = 2 — 6C, provided < ^ < 1/3, we obtain from the energy 
balance equation (27): 

d 112 , 211 ||2 , A 6C\ ^ l^^^l^+C^l/^-(:/2n II1+6C 

-(l|w||o + a ||w|L+A|||wr) < ^^3^. 11^11/ • 



This implies 



(29) 



dt 



where it is denoted 



D = 



/3C5/6+cC^l/3-C/2(^«)"'V4 




(30) 



e<((el*=o)'/'-'^ + (l/2-3C)i^t) 



(1/2-3C) 



-1 



<^(^;C), 



(31) 



and for C = 1/6 an exponential one 



e<ek=oe 



Dt _ 



^(t;l/6). 



(32) 



For ( > 1/6, (29) yields bounds for ^, which are finite-time and thus not of interest. In 
view of (24), (31) and (32), the solution to equation (22) satisfies the inequalities 



|||w||| < (<^/A)^/(6^) ^ 
|v(x,t)|^(^(t.^)/^)i-i/(3o,i/2 < M^,t)\i3i^^^,t)i-e,i/2 = ||w||i/2 < |||w||| < {ip{t;C)/A) 



Note that the first index of the norm in the l.h.s. of this inequality is strictly positive for 

any t > 0. 

We have proved the following statement: 

Theorem 2. Suppose < C < 1/6 and s > 5/6 + C- Let the norm | v*^"^) |o-_s_,_i/2 
of initial condition v^^'^) of a solution to the Voigt-type regularisation of the force-free 
three-dimensional Euler equation (22) be finite for some a > 0. Then the bound 



holds true for the solution at any time t > 0. Here A is defined by (14), e = 2 ~ 6( 
and /3 satisfies inequality (26); (/? is defined by (31) and (32), where ^|(=o is obtained by 
application of formula (30) to w(x, 0), which is the result of application of transformation 
(9) to the non-truncated initial condition v^™^ 

To the best of our knowledge, the regularised equation (22) with the fractional power 
of the Laplacian has not yet been considered in Sobolev spaces in literature. It is easy 
to show incrementally that for s > 5/6 the solution of (22) belongs to Sobolev spaces of 
arbitrary high indices (provided the initial data and the forcing f are sufficiently regular). 
Such a proof can be started by multiplication of (22) by v and deriving a bound for || vHs/g- 
Scalar multiplying now (22) by (— V^)^/^~'^v, using the inequality |v|9/2 < C5/6||v||5/6 to 
bound the integral arising from the nonlinear term, it easy to show that ||v||5/3 is bounded 
at any t > 0. Higher-index Sobolev norms can be subsequently bounded similarly. 



|v(x,t)|^(^(t.^)/^)i-i/(3c),i/2 < (v?(t;C)/^) 



1/(60 
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5 Voigt-type regular isat ion 

of the Navier— Stokes equation 



In this section we prove global (in time) boundedness of Gevrey-Sobolev norms of three- 
dimensional solenoidal space-periodic zero-mean solutions (3) to the Voigt-type regulari- 
sation of the Navier-Stokes equation 

«'^(- V + 1^ + (v . V)v - = f + Vp, (33) 

where s > 0. The Voigt regularisation (s = 1) was previously investigated in [25, 36]; 
boundedness of Gcvrey-Sobolev norms of solutions for s = 1 was demonstrated in [16] by a 
method, different from the one used here. We show here that a milder regularisation, than 
the one considered in the previous section for the Euler equation, suffices to guarantee 
global regularity of solutions. In the first two subsections we derive the bounds separately 
for s > 1/2, where the smoothing term is excessively strong, and for the limit value s = 1/2 
— constructions in the two cases arc somewhat different. In the last subsection we show 
that the method [11] can be used for s < 1/2 to establish instantaneous development of 
analyticity of solutions for mildly regular initial conditions. 



5.1 Excessive damping: s > 1/2 

Assuming | v*^™^ |o.,s+i/2 < oo for a positive a, we make transformation (24) for 1 < e < 2, 
where the norm ||| • |||, equivalent to || • ||s+i/2, is defined by (25); P is supposed to satisfy 
inequality (26), implying |||w(x, 0)||| < oo. We need to derive a bound for |||w(x, t)|||. 
Substituting (3) and (24), we transform (33) (for f = 0) into the system of equations 

(1 + aVP*) I ^T-^ + lllwlll"^"'' w„ -^|||w||| 1 + M'^Wn 
\ at at J 

= -i ^(Wk • (n - k)) PnWn_k exp (/3\\\w\\\-' (|n| - |k| - |n - k|)) . 

k 

The real part of the sum of these equations scalar multiplied by Wn reduces to 

^ {\\w\\l + a^w\\l + A\\\w\f-') + 2iy\\w\\l 
= Im^(wk • (n - k))(w„_k • ^) (exp (^]]]w|||-^ (]n| - |k| - ]n - k|)) 

k,n 

-exp(^|||w|||-^(|n-k|-|k|-|n|))) (34) 

(this identity differs from (27) only by the presence of the diffusion-related term 2i/||w||^). 

Note that |e"'-e°"] < \a'-a"\^^^ for all a' < and a" < (since £ > 1). Consequently, 
the r.h.s. of (34) does not exceed 
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^ |wk||n - k||w„_k||w„| (2^|||w|||-^ ||n| - |n - k||)V- 

n,k 

- (27r)3|||w||| J^Joi-^)h/e{^)M^)d^ < (2^)3|||w||| '^°'6l/lAl3/(.+l/2)l/ll3/(2-.) 

(by Holder's inequality and definition (13) of the functions fq) 

(2/3)^/" 

(by the Sobolev embedding theorem, see (7)). We have assumed here 1/2 < s < 1, and 
we further demand 2 > e > 1/s. Then, by Holder's inequality and definition (13) of the 
functions fg, (35) is bounded by 

CiCi-,C._i/2(2/j)V- 

(27r)3|||w||| 

where we have denoted 

= CiCi_,C,_i/2(2/3)i/^(27r)-Va, k = min(l, {s - l/s)/{l - s)). 
We therefore obtain from (34), by Young's inequality. 



dt 

where 



2^2 V 4i/ 
Therefore, 



Ikllo + «lw||2 + ^|||w|f < (||w||^ + a^w\\l + ^|||w|f -^) e^^.^* = ^{t). (36) 
Since |||w||| < {Lp/Ay^^'^~^\ transformation (24) implies 

|v(x,t)|^(^(t)/^)-e/(2-e),i/2 < |v(x,t)|^|||w(x,t)|||-e,i/2 = ||w||i/2 < |||w||| < {ip{t) / A)^/ 

The first index of the Gevrey-Sobolev norm in the l.h.s. of this inequality is strictly 

positive for any t > 0. 

Thus, we have demonstrated 

Theorem 3. Suppose 2 > e > 1/s > 1. Let the norm | v*^"') |o- s_|_i/2 of initial condition 
v^"^) of a solution to the Voigt-type regularisation of the force-free three-dimensional 
Navier-Stokes equation (33) be finite for a cr > 0. Then the bound 



l/{2-e) 



holds true for the solution at any time t > 0. Here A is defined by (14), /3 satisfies 
inequality (26) and ip is defined by (36), where w(x, 0) is obtained from the non-truncated 
initial condition v*^™^ by transformation (9). 
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5.2 The critical damping: s = 1/2 

In this subsection we focus on the hmit value s — 1/2 and show that in this case bounds 
for the Gevrey-Sobolev norms can be obtained in a similar way. Instead of (24), we now 
make a substitution 



v„(t) = w„(t)exp (-/3|n|(l + |||w(x, t)|||)-^J . (37) 
III w(x, 0)111 is bounded uniformly over truncations of v*^'") provided 

0</3<a(l + ^|v(-)|^_^/2 + «'|v('°)|^,i)' (38) 

(under this condition a transformation of the non-truncated initial condition v^™) is also 
well-defined). 

As in the previous subsection, wc derive from the Voigt-type regularised Navicr-Stokes 
equation (33) (for f = 0) a system of equations governing the evolution of Fourier coef- 
ficients Wn(i) of w, and consider the real part of the sum of these equations, scalar 
multiphed by w„. The energy balance equation, analogous to (34), takes now the form 

^(llwllS + alwll?,,) + AfS ^JWt_||||w||| + 2.||w||? 
= Im ^(wt . (n - k))(w„_t . *;) (exp (/3(|n| - |k| - |n - k|)/(l + |||w|||)2) 



k,n 



- exp (/3(|n - k| - |k| - |n|)/(l + |||w|||)2)) . (39) 

Using the inequality |e°' — e""| < jo;' — q;"|^/^ which holds for any a' < and a" < 0, we 
obtain from (39) the energy inequality 



i (llwll? + aVilf,, + 4^ (ln(l + lllwlll) + |^)) + 2Hlw|i; 



< y2^X:|wk||n-k||w„_t||w.||k|''^(l+ |||w|||)-' < -^a^aT^Tr^ 



n.k 



< ^^ff^ llwlli (40) 



(27r)3a2 
Consequently, if 

(in addition to condition (38)), then for 

i{t) = + alw||?/2 + 4/31n(l + |||w|||) (41) 
we find, integrating (40) in time: 

rt*Wrtm^i«/' 4|||w(x.t)|||+3 , 4|||w(x.0)|||+3 , 

m < m + ^-2(i + |||^(,.,)|||). + 2(i + iiiw(x,o)iiir^j ^ + 

12 



Each of the three terms, constituting ^{t), is positive; hence we deduce from this inequahty 



1 + |||w(x,t)||| < exp(e(0)/(4/3) + 3/2), ||w(x,^)||o < ^/m + 6p. 
In view of (37), the two inequahties imply 

|v(x,i)|^exp(-^(o)/(2/3)-3),o < | v(x, t) |^(i+|||^(^_t)|||)-2,o < ||w(x,i)||o < ^J^{0) + 6/3. 

Thus, in the case of critical damping the following theorem holds: 

Theorem 4- Suppose a > and the norm | v^™) of initial condition v^"^^ of a solution 
to the Voigt-typc rcgularisation of the force-free three-dimensional Navier-Stokes equation 
(33) for s = 1/2 is finite. Then the bound 

|v(x,t)|^exp(-5(0)/(2/3)-3),0 < ^^(0) + 6/3 (42) 

holds true for the solution at any time t > 0. Here (3 satisfies inequality (38) and ,^(0) is 
determined by application of formula (41) to w(x, 0), obtained from the non-truncated 
initial condition v*^™^ by transformation (37). 

What happens if the initial data v^''^^ is non-analytic and thus Theorem 4 is inappli- 
cable? Like in the case of the Voigt-type rcgularisation of the three-dimensional Euler 
equation, it is possible to develop the theory of solutions to equation (33) in the Sobolev 
spaces. In particular, one can show incrementally that for s > 1/2 the sohition of (33) 
belongs to Sobolev spaces of arbitrary high indices (when the initial data and the forcing 
f are sufficiently regular). Since this question is not in the scope of our paper, we only 
present a brief sketch of derivation of the bounds. Multiplication of (33) by v demon- 
strates boundedness of Jq || v(x, r) Hfdr for any t > 0. Scalar multiplication of the equation 
by (— V^)v and the use of the inequality 

|(v • V)v • (-V')v| < IvIelVvlal W|2 < (C||v||?||v||?+, + i/||v||^)/2 

yields 

|(«^l|v||?+. + ||v||D<l|f||L. + (i + q|v||Dllv||^^^^ 

whereby 

o^'MUs + Ml < («lv(x, 0) + ||v(x, 0) II?) exp (a-' j\l + C||v(x, r) \\l)dT^ 

+ ||f(x,r)||L.exp + C||v(x, r')||D^r') dr. 

Bounds for higher-index Sobolev norms can be subsequently derived by a similar proce- 
dure. 
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5.3 Instantaneous development of analyticity for s < 1/2 

The remark, made in [16], that the Navier-Stokes-Voigt equation (s = 1) exhibits damped 
hyperbohcity, remains vahd for the milder regularisation for s > 1/2 considered in subsec- 
tion 5.1. In particular, solutions for s > 1/2 apparently cannot instantaneously acquire 
analyticity — at least, the method [11], revealing that solutions to the Navier-Stokcs 
equation are capable of this, is not directly applicable to equation (33). To see this, let 
us transform a solution to (33), following [11], using the relations 

v(x, t)^Yl w„(^)e-^*l"l+-•^ w(x, t)^Yl Wn(i) e'"•^ (43) 

riT^O n 

where /3 > is a constant. The evolution of the transformed Fourier coefficients satisfies 
the equation 



:i + a>|^^) - /3|n|wJ + Hn^w: 



n 



= -i^(wk-(n-k))P„w„_kexp(;gt(|n| - |k| - |n-k|)). 

k 

Scalar multiplying it by |np'''Wn for 7 > 0, summing over n and taking the real part we 
find 

= Im^(wk ■ (n - k))(w„_k ■ ^)|np^exp {/3t{\n\ - |k| - |n - k|)) . (44) 

k,n 

Thus, the strength of the regularising term for s > 1/2 gives rise to a problem: viscous 
dissipation is too weak to control the term /5Q;^||w||^y2+7+s) appearing in the l.h.s. of 
the energy balance equation (44). Still, for s > 1/2 finite-time Gevrey class Gi/^2—2s) 
regularity emerges instantaneously, and this can be established by the method [11]. 

To show that the method [11] works for < s < 1/2, we choose three quantities rji > 
such that rji + r]2 + rjs — u and 1/2 < 7 < 1. If s = 1/2, we also demand 

/3 < ml a" (45) 
in transformation (43). By Holder's and Young's inequalities, 

/3||w||?/2+, < r^illwll^^^ + (4r7i)-Vlw||^, (46) 
^«lw||?/2+7+. < ?72||w||?+^ + Qiq;^||w||^+^, (47) 

where 

r 0, if^<w«', 

— -, ^ -, r , otherwise. 

[ 2(1 -s) \2r]2{l-s)) 

Using Holder's inequality, the Sobolcv embedding theorem (7) and Young's inequality, we 
find that the r.h.s. of (44) is bounded by the sum 
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^|wk||n- k||Wn-k||Wn||n| 



27 



k,n 



< (27r) Vl6/(3-27) 



E 



Wnllnl e'"''' 



^|w„||n|2^e-'"-'' 



6/(1+27) 

l/2+7||wll^/2-'>' 



< (27r)-^C^||w||^Ci/2||w||3/2Ci-^||w||i+^ < (27r)-^C^Ci/2Ci_^||w||;/^+^||w||^;';- 

< (g2/2)||w||f + m\MU„ (48) 

where we have denoted 

Q2 = 2(27-i^^^"^^^ rC',Cv2Ci_, 



4(27r)3 

Relations (44)-(48) imply the energy-type inequaUty 

d 
dt 

Therefore, 

df 



^(||w||^ + «^||w||^+J < {2m)-'p'\\wr^ + 2gialw||^+, + Q^llwllf +^^)/(^^-^). 



/ ..^ , ^ ^(1+27)7(27-1) 



dt 
where 

^{t) = ||w||^ + alw||2/2+^, ? = max((2r^i)"'/S', SQi)- (49) 
Integrating this inequality, we obtain a bound 

m < ((e(0))-(^-^/^)"^ - (g2/g)(e^*(^-^/^)-^ - 1))"^""'^'^ ^ ^(t), (50) 

valid for 

'<'-"^"'(^ + g.K(o)^(7-v2' )- (51) 

From transformation (43), relation (49) and this bound we infer 

Theorem 5. Suppose 0<s<l/2<7<l and the initial condition of a solution to 
the Voigt-type regularisation of the force-free three-dimensional Navier-Stokes equation 
(33) belongs to the Sobolev space if^+i/2(T^). Then for t < t^^ the solution satisfies the 
bound 



|v(x,t)|^t,^ < ^Jifit). 

Here a positive constant (3 satisfies inequality (45), if s = 1/2, and is arbitrary otherwise, 
if and t^: are defined by formulae (50) and (51), respectively, where ^(0) is determined 
applying (49) to w(x, 0) v^™). 

This shows that analyticity of solutions to the Voigt-type regularisation of the Navier- 
Stokes equation (33) for s = 1/2 emerges instantaneously, provided the initial conditions 
are in the Sobolev space if^+i/2(T^). For small values of the regularising parameter a the 
bounds ip{t) and are uniform in a (note that q is independent of sufficiently small a). 
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6 Concluding remarks 



We have explored a new approach to derivation of inequahties for Gevrey-Sobolev norms 
of solutions to evolutionary partial differential equations, in which a suitable nonlinear 
transformation of a solution in the Fourier space introduces a feedback between the norm 
of the transformed solution and the first index of the norm (6). Several examples of 
application of our technique were discussed. 

We have proved that if initially a three-dimensional flow is analytic, then analyticity 
in spatial variables is preserved by the Euler equation on the interval [0, t*) (16) — this 
is imphed by the bound (19) for Gevrey-Sobolev norms of a solution. (Alternatively, 
a bound for a Gevrcy-Sobolcv norm of the solution can be obtained [34] from finite- 
time bounds for solutions to the Voigt regularisation, which are uniform in the small 
parameter in the regularising term (a in (22)); such uniform bounds can be derived [34] 
by application of the method similar to the one used in [21].) 

Pivotal to our technique is transformation (9), which results in emergence of a new 
term, 

/3£ I n| 1 1 W (X, t ) 1 1 W„ — 1 1 W (X, i) 1 1 ,+3/2 , 

in the equations (10) and (21) governing the evolution of a solution (more precisely, of its 
transformed Fourier coefficients) to the Euler and Burgers equations. This term represents 
a new pseudodifferential mildly diffusive operator, controlling the norms in Hs+3/2', it is 
analogous to the term 

-f(t)||(-v2)'-+^/V||2 

appearing in literature on analyticity of solutions to PDE's of hydrodynamic type [17, 22]. 
(Our assumption s < 1/2 is technical, similar bounds can be obtained in other Sobolcv 
spaces for larger s in a slightly different way.) The order of this operator cannot be 
increased without acquiring difficulties in deriving bounds for the exponent in the energy 
balance equation (11). Evidently, by the same construction one can derive an identical 
bound on the interval [0,t*) for solutions to the Navier-Stokcs equation. This bound is 
uniform in viscosity and can be employed for a study of convergence of solutions to the 
Navier-Stokes equation to the solution to the Euler equation when viscosity tends to zero 
(such a study is beyond the scope of the present paper) . 

Solutions to the Burgers equation are known to develop shocks at finite time [14] 
due to intersection of characteristics; this rules out global (in time) analyticity in spatial 
variables of its solutions. That our bounds for the Gevrey-Sobolev norms of solutions to 
the Euler and Burgers equations are identical suggests, that the bounds are rough. Other 
nonlinear transformations of solutions can probably yield more accurate bounds. A clear 
drawback of our technique is associated with its relative simplicity: we cannot use it to 
demonstrate persistence of analyticity of solutions to the Euler equation, i.e. the fact that 
the solution is analytic on any interval < t < tcj where it is continuously differentiable 
[2] (persistence of non-analytic Gevrey-class regularity of space-periodic solutions to the 
Euler equation is shown in [17]). 

We have also investigated the Voigt-type regularisations (22) of the Euler equation, 
involving the regularising term a'^d/dt {—V'^Yv. and proved boundedness of Gevrey- 
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Sobolev norms of solutions to the regularised Euler equation for s > 5/6. Thus, a regu- 
larisation, milder than the one investigated in [21] (for s = 1), suffices to guarantee global 
regularity and analyticity of the solutions. Our technique enables us to decrease the order 
of nonlinearity of the energy-type inequality for the norm of the transformed solution; as 
a result, our bounds are global in time and exhibit a polynomial or, at most, exponential 
growth in time. 

Finally, we have explored the Voigt-type regularisation (33) of the Navier-Stokes equa- 
tion with the regularising term for s > 1/2, which is milder than the regularisation con- 
sidered in [16] (s = 1) and the regularisation of the Euler equation (s > 5/6) studied 
here. We have obtained global exponential in time bounds for Gevrey-Sobolev norms. 
For s — 1/2, both the method [11] and our technique can be apphed. Consequently, if 
the initial velocity is in the Sobolev space H^^i/2{T^), then the solution has a bounded 
Gevrey-Sobolev norm (50) on the time interval [0, t*) defined by (51). Furthermore, start- 
ing at any point tg, < to < t^, we obtain a global (in time) bound (42) (performing 
optimisation of the bounds in and other parameters is desirable). Unlike all other 
bounds for Gevrey-Sobolev bounds that we have derived in this paper, the bound for 
solutions to the Voigt-type regularisation (33) of the Navier-Stokes equation for s = 1/2 
does not deteriorate in time — neither the indices of the norm, nor the r.h.s. depend on 
time in inequality (42). 
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